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The dynamics of a single spin embedded in the tunnel junction (quantum point contact) between 
ferromagnets is addressed. Using the Keldysh technique, we derive a quantum Langevin equation. 
As a consequence of the spin-polarization in the leads, the spin displays a rich and unusual dynamics. 
Parallel configured and equally strong magnetic moments in the leads yield an ordinary spin pre- 
cession with a Larmor frequency given by the effective magnetic field. Unequal and/or non-parallel 
configured magnetization, however, causes nutation of the spin in addition to the precession. Our 
predictions may be directly tested for macroscopic spin clusters. 
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The interest in a number of techniques that allow one 
to detect and manipulate a single spin in the solid state 
remains tremendous both experimentally [TJ |2j [3j 0] and 
theoretically EJ EE QSJ E]. Being a crucial ele- 
ment in spintronics and spin-based quantum information 
processing, such studies are also of fundamental impor- 
tance. So far, most of efforts |3 HI GO 113 EI] have been 
focused on understanding the mechanism for the tun- 
neling current modulation, which is the hallmark of a 
single spin detection by using the scanning tunneling mi- 
croscopy (STM) [1, 2, 3 . Recently, the coupling between 
a single spin and supercurrent in Josephson junctions has 
also been studied [12 [T5] . 

Especially, a single spin nutation induced by an ac su- 
percurrent in a dc biased Josephson junction was shown 
for the first time in Ref. [12]. From the view point of 
single spin manipulation, such kind of nutation provides 
a significant implication for the control of spin dynamics 
electrically. The question is whether the spin dynamics 
(or even switching) can be realized in other types (i.e., 
non-superconducting) of leads, which are easily accessible 
experimentally. Recently we have examined this issue in 
normal conducting leads [14 . It was found that the spin- 
flip process of tunneling electrons is important to manip- 
ulate the spin and the flip-rate determines the efficiency 
of spin manipulation. Studies of local spin-dynamics in 
quantum dots between ferromagnetic leads were recently 
reported [131 US], however, these studies were concerned 
with the time-dependent effects and the noise of a local 
spin under the influence of stationary external fields. 

In this paper, we study the spin dynamics of a sin- 
gle spin embedded in a tunneling junction between two 
ferromagnetic leads. A quantum Langevin equation is 
derived for the single spin dynamics. Through the re- 
sulting equation we show that the tunneling between the 
ferromagnetic leads converts the electric field, e.g. bias 
voltage, into an effective magnetic field. The resulting 
effective magnetic field, however, depends on the relative 
orientation between the magnetization in the two leads. 
Parallel (i.e., FM-type) alignment of equally strong mag- 



netization in the two leads yields a shift of the Lar- 
mor frequency, which scales as the squared ratio between 
the induced and the external magnetic fields. The FM- 
type alignment but with unequal magnetization or anti- 
parallel (AFM-type) alignment in the two leads, on the 
other hand, leads to nutations in the precession of the 
spin about the external field. 

The model system we consider consists of two ferro- 
magnetic leads coupled to each other by a single spin S. 
We assume that the tunnel junction is formed by a quan- 
tum point contact between the leads, so that the mag- 
netic fields generated at the tips of the magnetic leads 
can be neglected. We also neglect the direct interaction 
of the spin with the two leads. Then, the system Hamil- 
tonian can be written as 

H = Hl + H r + H s + Ht. (1) 

The first two terms H l /r = J2kaeL/R, £ k<jc\ a c k< j describe 
the electrons in the leads, where an electron is created 
(annihilated) in the left/right (L/R) lead at the energy 
Efccr by c' ka (cfe CT ). Henceforth, we assign subscripts p (q) 
to electrons in the left (right) lead. The Hamiltonian for 
a free spin S in the presence of a magnetic field B is given 

by 

H s = -gii B B ■ S, (2) 

where g and [Ib are the gyromagnetic ratio and Bohr 
magneton, respectively. The two leads are weakly cou- 
pled via the tunneling Hamiltonian 

7i T ( C lai T a S a/3 + TiS ■ CT a [i]c q f3 + H.C.) . (3) 

pqaft 

Here, cr a p is the vector of Pauli spin matrices, with spin 
indices a, j3, whereas T a is the direct tunneling rate be- 
tween the spin-polarized leads, and T\ is the tunneling 
rate between the leads modulated by the local spin. In 
this model we neglect spin-flip transitions in the direct 
tunneling between the leads. For convenience we take 
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the respective amplitudes to be momentum independent 
(although it is not required) . Typically, from the expan- 
sion of the work function for tunneling T\/T a ~ J/U [7], 
where J is the exchange coupling strength between the 
transport electrons and the local spin and U is a spin- 
independent tunneling barrier. We further allow a weak 
external magnetic field B ~ 10 2 — 10 4 Gauss, which is ap- 
plied along the z-direction in the z-x plane perpendicular 
to the electron tunneling direction (y axis). 

When a time-dependent voltage bias is applied across 
the tunneling barrier, such that V(t) = Vd c +V ac cos(ujot), 
where Vdc and V ac are the dc and ac components, and ojq 
is the frequency of the ac field, a dipole will be formed 
around the barrier region through the accumulation or 
depletion of electron charge. This process results in the 
time dependence of single-particle energies: E p = e p + 
Whit) and E q = e q + Wii(t), with the constraint Wl(£) — 
Wu(t) = eV ac cos(oJot). However, the occupation of each 
state in the respective contact remains unchanged and 
is determined by the distribution established before the 
time dependence is turned on. Therefore, the chemical 
potentials on the left hl and on the right lead fiR differs 
by the dc component of the applied voltage bias, /jl — 
fiR = eVdc- The tunneling junction with the spin then 
has two time scales: The Larmor precession frequency of 
the spin wl — 9^bB and the characteristic frequency uiq 
of the ac field. 

The model is gauge transformed by 

fj = e - 1 ti [HL+W L (t')]N L dt> e -i tf o [» R +W R (t')]N R dt' ^ 

with Nl/r = J2k<j£L/R c l:cr c k<Tj and we thus obtain the 
model K — Kl + Kr + K$ + Kt, where K L i R — 

J^kaeL/R^kacl^Cka, £fc<x = Efecr - (J>L/R, K S = Hs, and 

K T = E (iXpc q 0e im + H.c.) . (5) 

pqa{3 

Here, we have introduced the notation T a p(t) = T a 5 a p + 
TiS(t) • cr a0 and <j>{t) = e J t * [V dc + V ac cos cj t']dt'. 

We now derive the effective action via the Keldysh 
technique [T7j. If all external fields are the same 
on both forward and backward branches of the 
Keldysh contour (C), then the partition function Z = 
tr Tq exp [— i § c Kx{t)dt] — 1, where the trace runs over 
both electron and spin degrees of freedom. We take the 
partial trace in Z over the lead electrons (bath) to ob- 
tain an effective spin action. In the present situation, 
this action represents the interaction of the magnetic 
spin with a non-equilibrium environment. The tunnel- 
ing contribution to the resulting spin action reads i8S — 
-(1/2) § c § c (T c K T (S(t),t)K T (S(t),t'))dtdt', much in 
the spirit of Refs. [18] . 

For brevity we put A a p = J2 pq c P a c qP- The tunneling 



Hamiltonian of a voltage biased junction then reads 



K T (S(t),t) = J2(T a p(t)A aP e 1 ^ + H.c). (6) 

a/3 



For magnetic leads, the correction to the effective action 
for the spin dynamics is thus given by 



iSS = -i I I fap^Vapit^t^fp^t'y^-^'^dtdt 1 ,(7) 
Jc Jc 



where V a/3 (t,t') = -i(T c A Q/3 (i)v4 T Q/3 (t')>- 

Performing standard Keldysh manipulations, defining 
upper and lower spin fields S u /' residing on the for- 
ward/backward contours and reducing the time ordered 
integral over Keldysh contour to the integral over forward 
running time at the cost of making the Green functions 
(GF) V a p a 2 x 2 matrix, for each spin combination a/3. 
We then perform a rotation to the classical and quantum 
components 



S c = (S u + S')A S 9 = S" - S l , S c • S 9 = 0, (8) 



which makes the matrix Green function uniquely de- 
termined in terms of the retarded/advanced component 
K'pW) = T#(±t T t'){[A a0 (t),Ai (t')]), and the 
Keldysh component T>^(t,t') = -i({A af3 (t), A^t 1 )}). 
The procedure leads to SS = 8S C + SS q , where 



5S C = J [SI(t)KP(t,t') + Sf(t)Kg\t,t')S^(t')}dtdt' 

(9) 

and 



ss q = Js?(t)K§\tj)s]V)dtde, (io) 



where summation over repeated indices i,j — x,y,z, is 
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understood. Here the kernels are given by 

a 

+V^(t',t)e- lWt) -' t ' it ' )] ] (11) 
= -2T x e{t - if) ° Z Mf(Z p(T ) - f(£ ga )] 



pqa 



x sin [(£ pff - i q «){t - t') + 4>{t) ~ </>(*')] 



a/3 



-i[<P(t)-<l>(t')h a 3_ 



+Vl a {t',t)e r 0a 

= nfo(t - f ) J2 < p {[.f{t pa ) - f&p)] 

pqaft 

xe i{ipc,-iaf}){t-t')+i[<t>{t)-<l>{t')] 
xe - i («p/J-« g =)(t-*')-i[0(O-0(*')]} cr J 



(12) 



(13) 



a/3 



pqa ft 



-2/fea)/(e 9 /3)]e 4( «--^ )(t -*' )+l ^ (t) -^' )] . 

Here /(£) is the Fermi distribution function /(£) = 
l/[exp(£/fcsT) + 1]. It follows that the only non- 



K {2) K {2) 



vanishing components of are K v \} 

and Kyx — Kxy • 

To properly describe the dynamics of the spin, we em- 
ploy the path integral representation for the spin fields. 
So in addition to the terms, — § Hs(t)dt, the action for 
a free spin also contains a Wess-Zumino-Witten-Novikov 
(WZWN), Swzwn, which describes the Berry phase ac- 
cumulated by the spin as a result of motion of the spin on 
the sphere. We generalize this action for nonequlibrium 
dynamics within the Keldysh contour formalism, which 
can be expressed as [12] 



i j dts 9 -(s c xd t s c ). 



Swzwn = ^ / dt S 9 • (S c x d t S c ). (14) 
The total effective spin action is given by: 
S c ff = Swzwn + g^B J dtB ■ S q (t) + 5S C + SS q . (15) 



the quadratic term in SS q . The total effective action is 
rewritten as: 

S cS = Swzwn + g^B J [B + £(*)] • S q (t)dt (16) 

+ J [smK^Kt^+sumifit^s^jdtdt', 

where the fluctuating random magnetic fields satisfy the 
correlation functions 



(17) 



Equation ( |16| ) constitutes the central formula for the 

(2) 

following analysis. Notice that the kernel K\- , which 
connects the quantum and classical spin fields (see 
Eq. (16l), represents the effects of electron degrees of 
freedom. It has the following approximate behavior in 
time: 



where 77 is an infinitesimal and fi is the energy scale at 
the order of voltage bias. A similar type of approximate 
behavior has also been observed in the case of a me- 
chanical oscillator coupled to two Luttinger liquids [19] . 



(2) 

Equation (18) suggests that the kernel Kh is peaked at 



t — t' = wmlc oscillates at larger times. The character- 
istic time is set by l//z, which also measures the width 

of the peak around t — t' = 0. Because of the oscillating 

(2) 

nature of at the large time scale, the corresponding 



time integral on the right-hand side of Eq. ( 16 1 is domi- 
nated from the time range of 1/fi around t—t = 0. When 
1/fj, <C 1/Vl (i.e., /1 3> Wi), we are in the regime where 
the spin dynamic processes are much slower than those 
of the conduction electrons. Since the energy associated 
with the spin dynamics, Hujl ~ 1/ieV, while the typical 
energy for the electronic degrees of freedom at the order 
of 1 meV is routine, the above regime is easily accessible. 
Under this condition, it is reasonable for us to use the 
approximation S c (t') » S c (t) + (f - t)dS c (t)/dt. The 
variational equations dS c s/SS g (t) = then yield 

^ = a(t) ^ x n + g t i B n X [B cff (i) + £(*)] , (19) 

where we have put S c (i) = Sn(t), whereas a(t) = 
S J K< 2 )(i, t')(t-t')dt' , with the 3 x 3 matrix K (2 \t, t') = 



{Kg\t,t')}i, j=x , y , z , and Beff(t) = B(f) + B&(t) 



( 2) 

As seen from Eqs. ( 11 1-( 13 ), the first three terms on B\ n ' d (t). We expect that the Langevin term (£(f)) is 

suppressed at frequencies much lower than the exchange 
interaction field in the leads. 

First we note that a(t) ~ STiUIq/D at zero tempera- 
ture, where 2D is the band-width in the leads. Hence, for 
large D, which is reasonable for metals, a is negligibly 



the right hand side of Eq. ( 15 1 are real, which deter- 
mine the quasi-classical equation of motion, while SS q is 
imaginary, which stands for the fluctuations of the spin 
field S q . This means that the quantum effects have in- 
deed been included even in the semi-classical approxima- 
tion. We perform the Hubbard- Stratonovich transforma- 
tion with an auxiliary stochastic field £(i) to decouple 



small such that first term to the right in Eq. (19) drops 
out. It is reasonable to believe that this would be true 
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also for finite temperatures. This result is different from 
the case of a dc biased superconducting tunnel junction, 
where a(t) is finite leading to spin nutation [12] . 

The magnetic leads induce the field B^H = B^lz with 

m B B^ d (t) = J KP(t,t')dt' 

a nm 

X J m (eV ac /uo)(eVdc + muj ) skiLJ Q (n - m)t, (20) 

where Nhalita is the density of states of the spin a sub- 
band in the left/right lead, whereas J n (x) is the nth 
Bessel function. This induced magnetic field contributes 
to the spin dynamics whenever, at least, one of the leads 
is spin-polarized. 

(2) 

The second induced magnetic field, B- n j(i) is defined 
such that n x gn B B$ d (t) = n x S fKW(t,t')dt'n, and 
we find that this field can be written as 

gf x B B^ d (t) = S j[K^(t,t')(n y k-n x y) 

-[Ki 2 J(t,t')-KW( t ,t')}n z z\dt' (21) 

Here, the longitudinal component is proportional to 

[K^(t,t')-Kf)(t,t')]dt' = 
= 2irTl(N L] - N Ll )(N m - N Rl ) £ J n (eV ac /tu ) 

urn 

xJ m (eV ac /ujo)(eV d c + muj ) sinw (n - m)t, (22) 

which thus contributes to the spin dynamics whenever 
both leads are spin-polarized. The transverse component 
of the induced field, which is proportional to 

K^(t,t')dt' = 

-- -2nT 2 (N L] N Rl - N Ll N m ) £ J n (eV ac /io ) 



the leads and describe the effect of the current flow on 
the local spin dynamics. The field B[^ arise whenever 
there is a spin imbalance in, at least, one of the lead. 
By its presence the local spin exerts a precession about 



its local direction. The z-contribution of the field B 



(2) 
ind' 

given in Eq. (22), is finite only when both leads are 
magnetic, and its effect on the local spin is analogous 



to that of Bp^ j. T he transverse contribution from B; nd . 
given in Eq. ( 23 1 , which is of main interest in this pa- 



(i) 



per, requires that the magnetic moments of the leads 
are non-parallel and gives maximal effect when they are 
anti-fcrromagnetically (AFM) aligned. Caused by the 
misaligned magnetic moments of the leads, the tunnel- 
ing electron has to undergo spin-flip transitions in order 
to tunnel between the leads. The spin-flipping tunnel- 
ing electrons tend to change the local magnetic field in 
the vicinity of the local spin, such that its local direction 
is slightly altered. In the case with AFM aligned leads, 
while the local spin tends to line up with the magnetic 
moment of the source lead, an ac component in the bias 
voltage tends to cause wobbling of the spin, as is shown 
below. 

The classical equation of motion should be consistent 
with the parametrization of the spin on the unit sphere, 
Sn = S(cos <^> sin #, sin <j> sin 6, cos 9) . Letting the external 
magnetic field B be oriented along the z-axis, we find the 
equations 



dt 



-g^ B [B(t)+Bgl(t)} 

+S [ [Ki 2 J(t,t') - Ki 2 J(t,t')]dt' cos9, 



(24) 



f t = S I KW(t,t')dt'smO. 



Defining A nm (V dc ,V ac ) such that S J (t,t')df = 
Enm A nm (Vdc, V ac ) cos w (n - m)t, and letting the spin 
initially at time to — be oriented at an angle 9q relative 
to B we explicitly find the polar coordinate 



x J m (eV ac / uo) (eV dc + muj ) cos v (n-m)t, (23) Q(t) = 2arctan 



tan 



vanishes for magnetic leads whenever N^Nm = 
NliNr^, that is, when the magnetizations of the leads 
are equal and in parallel configuration. Thus, under those 
conditions, the resulting motion of the spin will be in- 
fluenced by the magnetizations of the leads. More in- 
terestingly though, whenever the magnetizations of the 
leads are non-equal and/or in a non-parallel configuration 
(including anti-parallel) we expect the spin dynamics to 
become significantly modified, since then the equation 



o tt (a sinuj (n - m)t 
— I I exp^ A nm — 



Generically, when a single spin is subjected to a uni- 
form magnetic field, the spin azimuthally precesses with 
the Larmor frequency ujl- For the tunneling junction 
between ferromagnets, when a dc voltage bias is ap- 
plied the spin precesses around an effective magnetic field 
B ff — B + B^ 1 ) + B^ 2 ) with a shifted Larmor frequency 
u>l. When an ac voltage bias is applied, Eq. (24) shows 



of motion (19 1, contains also a non-trivial ^-component 



which is induced by the non- vanishing transverse compo- 
nent of the induced magnetic field. 

The contributions expressed in Eqs. ( 20 )-( 23 1 are 



that the spin also exhibits polar (9) modulations, when- 
ever the magnetizations of the leads arc non-equal and/or 
not in parallel configuration. In particular, when the ex- 
ternal magnetic field B(t) = 0, the effective field in Eq. 



(119} becomes B off (i) = B&(t) + B\^> d (t). Hence, we ex 



(2), 



caused by the spin imbalance, or spin-polarization, in pect that the dynamics of the spin in absence of external 
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gle variation equals that of the ac field, whereas the Lar- 
mor frequency of the spin precession scales as the squared 
ratio between the external and induced magnetic fields in 
a dc case. The induced magnetic field could be probed 
with superconducting quantum interference devices. 

We thank A. V. Balatsky and I. Martin for useful dis- 
cussions. This work was carried out under the auspices 
of the NNSA of the US DOE under Contract No. DE- 
AC52-06NA25396 and supported by the LDRD and BES 
programs at Los Alamos. 



FIG. 1: (Color online) The different qualitative behavior of 
the spin S c in presence of an applied ac voltage bias when the 
magnetizations of the leads are equal and FM-type (dashed), 
and AFM-type (solid) aligned, a) The polar displacements 
of the spin dynamics. Here, we have defined spin-polarized 
density of states [2D] N L(r/R<T = N (l + rj lT p L/R )/2 and r? T/i = 
±1, where TiJV ~ 0.1, p L = 0.9, p R = p L (FM-type) and 
p R = — pl (AFM-type), and ujo = 3, which gives the induced 
field g/is-Bind ~ O.OleVb- b) The resulting generic spin motion 
on the unit sphere. 



magnetic fields is analogous to that in the presence of ex- 
ternal magnetic fields. In Fig. [l^i) we show the dynamics 
of the z-componcnt of the spin, S^, when the magnetiza- 
tions of the leads are aligned in FM-type (dashed) and 
AFM-type (solid) configuration. In the plot, we only 
illustrate the first harmonic ldq, e.g. n — m = 1 and 
^ ! A nm constant, for a pure ac voltage, since the 
complete spin dynamics is expected to be very compli- 
cated. Clearly, the polar angle 8 exhibits a modulation 
with a dominant u>q harmonic. The resulting dynamics 
of the spin can to much extent be compared with that of 
a rotating rigid top, as schematically illustrated in Fig. 
[l]}) . Analogous to a classical spinning top, the spin wob- 
bles along the polar direction in addition to the azimuthal 
rotations. The spin-polarized current leads to a full non- 
planar gyroscopic motion (nutations) of the spin much 
like that generated by applied torques on a mechanical 
top, i.e. magnetically induced nutations. Similar dy- 
namics is expected to occur even for S = \ since in a 
spin coherent state the Schrodinger equation is essentially 
classical. The effect found in the present system is much 
similar to that of a single spin embedded in a Josephson 
junction |12j . although the origin is of different nature. 
In addition, the proposed system is experimentally more 
accessible because such challenges as extremely low tem- 
peratures and device engineering of a Josephson junction 
can be avoided here. 

In conclusion we find non-planar motion of a spin em- 
bedded in a tunnel junction between fcrromagnets in 
presence of an alternating current, whenever the fcr- 
romagnets are unequally strong and/or in non-parallel 
alignment. The characteristic frequency of the polar an- 
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